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r^ . Abstract. The hot star wind momentum problem rj = Mvca/(L/c) » 1 is revisited, and 

II ■ it is shown that the conventional belief, that it can be solved by a combination of clump- 

;, . ing of the wind and multiple scattering of photons, is not self-consistent for optically thick 

C/3 ■ clumps. Clumping does reduce the mass loss rate M, and hence the momentum supply, re- 

quired to generate a specified radio emission measure s, while multiple scattering increases 

T '~j ' the delivery of momentum from a specified stellar luminosity L. However, in the case of 

5_( ' thick clumps, when combined the two effects act in opposition rather than in unison since 



clumping reduces multiple scattering. From basic geometric considerations, it is shown that 
this reduction in momentum delivery by clumping more than offsets the reduction in mo- 
mentum required, for a specified e. Thus the ratio of momentum deliverable to momentum 
required is maximal for a smooth wind and the momentum problem remains for the thick 
clump case. In the case of thin clumps, all of the benefit of clumping in reducing r] lies in 
reducing M for a given e so that extremely small filling factors / x lO""* are needed. 

It is also shown that clumping affects the inference of M from radio e not only by 
changing the emission measure per unit mass but also by changing the radio optical depth 
unity radius i?iad, and hence the observed wind volume, at radio wavelengths. In fact, for 
free-free opacity oc n^, contrary to intuition, Ry^i increases with increasing dumpiness. 

Key words. Circumstellar matter ~ Stars: mass-loss - Stars: winds, outflows - Stars: Wolf- 
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1. Momentum problem 

If one infers the mass loss rate M for hot massive (especially Wolf-Rayet) stars from the radio 
emission measure e - ^ rP'dV, using a smooth spherical wind model, one finds that the wind 
'momentum' rate p - Muca = rjL/c involves 77 » 1, where L/c is the radiative momentum 
outflow rate (Cassinelli & Castor n"973> . Insofar as such winds are believed to be radiatively 
driven, this poses a 'momentum problem', the solution of which has long been a hot topic in 
the field (Barlow et al. nMTl Abbott et al. [T986l Cassinelli & van der Hucht fTWTl Willis [TWTl 
Lucy & Abbbot lD^ Springmann lD^ Springmann & Puls l 19951 Gayley, Owocki, & Cranmer 
119951 Owocki & Gavlev lT999> . Estimates of 77 vary according to assumptions (e.g., arguing for 
a high value of L) but values of t] ranging up to nearly 100 are mentioned (Hamann & Koesterke 
I1998> . There are two main strands of argument quite widely believed to combine to solve the 
momentum problem, one being mainly observationally driven and the other mainly theoretical. 
The values of M associated with these large 77 are those inferred from a smooth spherical 
wind density model, the radio emitting material filling the volume. The contribution to e from 
any volume element AV is Ae x rp-AV oc M'/AV. If, however, the material is clumpy, filling only 
a fraction / - {nf- l{rp-) of the volume, then Ae is enhanced by a factor 1// for a given M. The 
mass loss rate M required to generate an observed e thus scales as M oc /'^^ in clumpy winds. For 
strong clumping (/ <K 1), this ameliorates the momentum problem, though the / = 10 '^ required 
to reduce Mvo:, by a factor of 100 seems very unlikely, so this clumping efifect alone cannot be the 
complete answer (e.g., Nugis & Lamers 120001 cite clumping corrected mass-loss rates yielding 
?7 a; 6). For example, making clumps very small increases their radio optical thickness and may 
make optically thin emission measures irrelevant. There is extensive observational evidence for 
large scale clumping in WR winds: the presence of narrow emission features moving out on 
broad wind emission fines (e.g., Robert et al. 119891 Robert et al. ll991l Moffat & Robert [TWTl 
Kholtvgin ll995t : broad band photometric and polarimetric fluctuations (e.g.. Brown et al. 1 19951 
Li et al. l2000> : and the absence of strong electron scattering wings (which scales as (n)^ rather 



than (n^) (HilUerlMB)- 

On the theory side, it has long been recognised that the limit p - Muoc ~ L/c is only true 
if (all) photons are scattered once only. If the wind scattering optical depth is high, the photons 
can, loosely speaking, be scattered 'back and forth' across the wind delivering momentum of 
up to 2/iv/c at each scattering (for thick clumps) until v is progressively dissipated by Doppler 
reddening at each momentum-delivering scattering on the moving matter The nature of this mul- 
tiple scattering has been described with progressively greater insight over the years. In particular, 
Gayley et al. J1995> showed that scattering back and forth across the entire wind is not required. 
Instead, the momentum is delivered in a series of random semi-local scatterings of photons as 
they diffuse outward, provided successive scatterings involve long enough paths to sample dif- 
ferent matter velocities. The essential feature is that of the large scattering optical depth t, which 



Send offprint requests to: john@astro.gla.ac.uk(JCB); li@astro.gla.ac.uk(QL) 



Brown, et al.: Wind momentum problem 3 

enhances the momentum deUvery rate to tLjc (e.g., Friend & Castor [T^83l Kato & Iben ll9921 
Netzer & Ehtzur 0*9931 Gayley et al. I1995t . because the diffusive delivery scales with the num- 
ber of scatterings Ns as nJ , while N^ - t^- Since the predominant driver is via the large 
opacity/cross-section associated with lines, Gayley et al. J1995> : and Owocki & Gayley <1999> 
have suggested that the issue is not so much a momentum problem as an opacity problem. 

The massive WR winds are still believed to be driven by line opacity (e.g., Lucy & Abbott 
I1993> . Unlike the less massive winds of OB stars, the WR winds have significant ionization 
gradients, that can substantially alter the hne opacity distribution with radius in the flow (Herald 
et al. 120001 Vink, de Koter, & Lamers EOOO> . Consequently, as the photons move away from 
the star, interact with a certain line opacity that exists at some radius r in the flow, and then 
escape, the photons encounter a new line opacity distribution at a different radius. Consequently, 
if there are gaps in the line distribution at one radius, those gaps can be filled by a different line 
distribution that exists in another part of the wind flow. The opacity problem then represents how 
effectively all of these gaps are "filled". 

Photon escape at gaps in the hne frequency forest reduces the flux mean opacity (or flux mean 
cross section cr per particle in our formulation) used in the gray approximation. The maximum 
that can be achieved by multiple scattering is reached when the number of random scatterings 
is so great as to Doppler shift the photons down to near zero frequency, the maximum Doppler 

1 /2 

shift per scattering being of order v/c for wind speed v. This requires t ^ Ng ~ c/v (^ 100), 
implying Mvca ~ L/c x c/uco or Mi'^/2 ^ L/2 which is the energy conservation limit. Available 
calculations of multiple scattering with real opacities can yield rj gains of order 10, that may 
explain some WR winds, but not the more extreme cases in which rj ~ 10^ is required. 

Since reduction of M (for a given e) by clumping and increase of momentum delivery by 
multiple scattering can each offer a factor of order 10 reduction in the momentum problem, there 
seems to be growing widespread belief that the momentum problem can be laid to rest (e.g.. 



Conti ll995t . However, this involves the tacit assumption that these two factors can operate inde- 
pendently and constructively, the impact of clumping on the effectiveness of multiple scattering 
never having been addressed (Hillier & Miller [T999b (although Shaviv ( I1998t has discussed the 
related topic of how opticaUy thick clumps increase the Eddington luminosity). Here we show, 
using simple geometric arguments, that this assumption is incorrect in the case of optically thick 
clumps, and that clumping, while reducing M, also reduces r, so making multiple scattering less 
effective. Essentially this is because clumping reduces the number of scattering centres compared 
to scattering off of atoms and also reduces < n >, for a given s. (Note that when discussing the 
effects of clumping it is essential to keep in mind that the observed s is held fixed. This fact is 
sometimes overlooked). 

We find quantitatively that, for thick clumps, the reduction in multiple scattering momen- 
tum delivery more than offsets the reduction in momentum required, the nett effect being that 
clumping worsens the momentum problem rather than solving it. 
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2. Single clump 

To illustrate the point, we first consider one thick scattering clump of mass M composed of 
atoms/ions of mass m. This is taken to have very high internal optical depth in the line-driving 
wavelength range so the clump as a whole is the scattering centre. Since we are not concerned 
with the wind speed profile v{r) but only with the final wind speed and momentum, we here 
approximate clumps as moving radially with speed v ^ Va, and to have the shape of a conical 
slice of radial thickness 6 and solid angle Q., the volume of the cone being r^Q.5 at distance r. 
We assume the clump to be optically thin at radio wavelengths, so its radio flux depends on the 
emission measure, but optically thick to lines for the stellar radiation at short wavelengths that 
are responsible for driving the flow. 

The emission measure s\ - rP'V (which measures the radio emission rate) of a single clump 
is 



{MImf eio 



«l - -J7^ = -f (1) 



where eio - (M/m)^/(R^^^il6) the clump emission measure at r = xR^d = ^rad, for Rmd is 
the radius of the radio photosphere which may be hundreds of times larger than the optical 
photosphere radius. Note that, for a prescribed si, M oc y/Q. for any chosen r, 6. 
On the other hand, the available rate of delivery of momentum is 

- ^ " n^ 

Pi, avail - ~"i~- \^> 

c 4n 
where we ignore scale factors of order unity due to the effects of gravity and of the backward 
scattering angular distribution function. The rate of momentum delivery required is 

Pl,req ~ —T- (3) 

rjv 



It follows that, for a given £i, the effectiveness of momentum delivery to a single clump is 

P 1 ,avail 



^ 



Anmu^cs^ 



qV' 



(4) 
6 ' 



P 1 ,^6^ 

This decreases as we make Q. smaller - i.e., as we make the clump clumpier - because the 
momentum /^ire^ required oc M oc yO., but the momentum />i,ava,; available oc Q, and the decline 
of the latter with O is dominant for a single clump. That is, making the O of a single clump 
smaller does reduce pi^req for a given ei but reduces pi ^j,^,, even more. So shrinking one clump 
of a given si makes it harder to drive it to terminal speed of known value Uco- 

Compressing the clump radially does help (in this single thick clump case) since reducing 6 
reduces Mv for prescribed ei , Q. 

3. Multiple clumps 

We now have to consider the effect of multiple scattering in the case of a multiple clump wind, 
since multiple scattering cannot occur in the case of an individual discrete clump. In doing so we 
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take all clumps to be optically thick in the UV but thin in radio, identical in size and mass, and 
use the gray opacity approximation, the clumps being driven by a spectral mean 'continuum' ra- 
diation flux. We are of course well aware that in reality there will be a distribution of clump sizes 
and masses. However, if one can prove that for any specific clump parameters, clumping reduces 
the benefit of multiple scattering, then the same must be true of the sum over any distribution of 
clump parameters so long as they remain thick. Put another way, the arguments that clumping a 
wind increases its emission measure, that multiple scattering increase momentum delivery, and 
that clumping reduces multiple scattering all derive essentially from geometric arguments and 
have nothing to do with the details of opacity or of clump size distribution (other than being 
thick). 

Retention of the conical slice shape described above, taking Q and 6 independent of r, means 
that the clumps expand in 2-D (transversely) rather than in 3-D, which is reasonable for a highly 
supersonic wind. The constant Q, 6 assumption also means that, for constant v, clumps occupy 
the same fraction (constant filling factor /) of the volume at all r. For spherical (3-D) clump 
expansion, linear radial expansion {5 oc r) would result, for constant v, in radial merging of 
clumps, which corresponds to an r-dependent filling factor / with / — > 1 as clumps merge. 
Situations with non-constant filling factor / - f(r) have been discussed by Nugis, Crowther, 



& Willis ( fT998b : HiUier & Miller ( fT999t : and Ignace et al. J2003> . We assume clumps are, on 
average, emitted uniformly over the stellar surface at a rate C in clumps per second. Then the 
space density of clumps at r is 

where we again approximate v - constant = Doo and the radio emission measure of one clump is 
again given by Eq. Q. Using Eqs. Q and (|5j, the total emission measure can be written 

s = r Anr^drNs, = -^ (-)' i- = ^^s,„ = N,ffSu„ (6) 

The last form is interesting, showing that the total emission measure s is just the initial emission 
measure of one clump si„ at r = R^^i times an effective number of clumps Ne/f = CRrndlv, 
namely that located in the range R^.^^ < r < IR^-^d. 

The mass loss rate M and the momentum delivery rate preg required are 

M = CM (7) 

and 

i?r,g = CMV. (8) 

where we neglect scale factors of order unity as we did in Eq.(|3. By Eqs. (|6j and (|8|i, we get the 
momentum delivery rate required for a given total wind emission measure e as a function of M, 
D., namely 

2 2 ^^ 

Preq = m V RradS—- (9) 
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We want to compare this with the momentum delivery rate available from multiple scattering of 
stellar photons and we take this to be given by (cf. Section 1) 



Pavail 



-nV^ 



(10) 



where r is the mean (gray approximation) line scattering optical depth of the wind and A^, = t^ is 
the number of scatterings of an escaping photon, t is also the 'covering factor' or the total solid 
angle of all the clumps as seen from the star divided by An - see Appendix. 

The wind optical depth for starlight due to lines treated in the gray approximation is (for 
individually thick clumps) 



n 



Qdr 



(11) 



where n is the density of scatterers, Q is the cross section of scatterers, and /?uv is the UV 
photosphere associated with the dominant Une driving. 

We choose to split the range into two sectors, r < d and r > d, where d is the distance at 
which an individual clump becomes optically thin radially. At r < d the individual scatterer is a 
clump of area r^Q and thickness 6, while at r > c/, it is an ion of area cr (the actual value adopted 
for cr being some frequency average over lines). Thus the optical depth integral expands to 

d ,\ Mcr RyjY 



-f 



NiryOdr + 



f 

Jd 



M CRuY 

N(r) — crdr - — 

m Anv 



Q. 



Ri 



-1 + 



'«^uv ^ 



where d satisfies 
Mcr 1 



m Qd^ 
and thus 

d 



1 



Ma- 
mRlyj 



1/2 



1 



Qi/2' 
Then equation J12> becomes 

xI/2 



CR, 



T = 



Attv 



Mcr 



Q'/2_Q 



Anv \ m I 



2- 



Mcr ^ 



1/2 



Q. 



1/2 



(12) 



(13) 



(14) 



(15) 



Consider the second term in expression (I15> . The ratio -^ is the total area of all the atoms 
in a clump, and R^^^ is the total area of a clump at r = ^uv- Since our analysis deals with 
individually thick clumps, we require -^ » ^uv^' ^° ^^ ^^^ S,^^ ^ -^ ^uv and neglect the 
second term in Eq. ( I15> and write 
C lMcr\^l^ , 



T - 



2nv \ m I 
If we express C in terms of e by using Eq. (|6}, then Eq. M6\ becomes 



In 



I m y 

\m) 



:6^''\ 



(16) 



(17) 



Now we get the available momentum using Eqs. ( I10> and (I17> in terms of M, Q. for a given 
e, namely 



. _ L I R,.,,ct'I^ 

Pavail I rx 



Km) ■ 



:5Q3/2. 



(18) 
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Comparing Eqs. (|5) and (I18> . we find a dimensionless measure of the effectiveness of mo- 
mentum delivery, as the ratio of momentum available to momentum required, as a function of M, 
D., for a given s, namely 



xrj _ Pavail 



Inv^cR^ad 



/a-y/2 
\m) 



rad 

M 



(19) 



For a given star (L) and a wind of given speed v and composition (m, cr), the expression in the 
bracket [ ] of Eq. ( I19> is constant, so we can vary the value of *P by varying the clump parameter 
combination Q/M. 

The essential result is that *? increases with increasing Q/M, i.e., with increasing clump cross 
section per unit mass (which is different from the single clump case of Eq. Q). To minimise the 
momentum problem (maximise W) for a given mass M (and thickness S), Q. should be as large 
as possible while for a given Q the mass M should be as low as possible with, in both cases C 
varying according to Eq. (|6} to ensure the correct s. If we change (e.g., increase) 6, ^ does not 
change but C changes (falls) to maintain fixed s. Consequently, to maximise *? we must make 
the clump mass small, the clump angle large, and the clump thickness large with correspondingly 
small C, all of these corresponding to minimising clumping. 

It is also of interest to express 'I' in terms of the volume filling factor / = (n)^/(n^) which 
can be expressed (with Vc = single clump volume - Q.r^6 at r) as 

where /q = £: and fr = ^ are sohd angle and radial filling factors respectively. Alternatively, 
the volume filling factor can be expressed as 

■' 4n \MJ ^ ^ 

Comparing Eqs. J19l ) and J21> . we see that, under our thick clump assumption, Q/M oc /'/^ 
so that *? oc Z'^"*, so that for any values of {s,6), decreasing / decreases the effectiveness of 
momentum delivery. 

All of the above shows that, contrary to conventional 'wisdom', in the case of thick clumps, 
clumping does not help solve the momentum problem but actually makes it worse. 

The case of a smooth wind can be considered a limit of the clumpy case as the clumps blend. 
However, there are infinitely many clumpy cases that approach the smooth case as the clumps 
blend and it is easier to evaluate ^F = ^Fj, for the smooth case directly using no(r) - M/{4-nr^mv). 
Then with subscript 'o' denoting the smooth case, we get 

Po,req = Mv, (22) 

m-^ 4;ri)'^Ki-ad 

- = r^. a4) 
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and 

L L Mcr 

Po,avail - -T„ - (25) 

c c AnmuKuY 

from which we deduce that 



L 1 o- 

c Anu-RnY m 



-^o^-T-n^- (26) 



Exphcitly comparing the momentum dehvery effectiveness for the clumpy and smooth cases 
we have, by Eqs. ( I19l l and ( I26t 






^l.,\j(T\\'^^ t(2Ruv) 



'uv 



M 






2\ 



Mcr 
mil 



1/2 



^^"^ (27) 



d 

which is clearly <k 1 for clumps which are initially optically thick. Note also that ^ involves the 
ratio of Q.{2Ruy)^/M and cr/m, respectively, the cross-sections per unit mass of clumps and of 
atoms, while ^'o involves only cr/m. Clearly, the continuous limit corresponds to Q.{2R\jv)^/M — > 
cr/m at t/ — > 2/?uv (cf. Eq. ( I13t ). the 'clumps' become individual atoms, becoming thick at that 
point. Note that d — > R\j\f is essentially the limit where driving approaches the smooth wind 
limit, equally applicable to optically thin clumps. 

Carrying this line of inquiry further, it is helpful to see how Eq. J19t for *F approaches the 
smooth limit ^'g. We require that t = r,, from Eqs. ( II6I 1 and J24> . that e = £„ from Eqs. (|6} and 
( I23> . and finally that f - I- These conditions are met for 



V 47T Rad 



^N.fffjn^h (28) 



— = Q(2«uv)', (29) 

m 

and 

167TmvR^„, 

6 = -r—^ (30) 

Mcr 
which, on substitution in Eq. jl9> gives ^P = ^Po (Eq. J26» as required. To interpret Eq. j29t 
physically, note that Mcr/m is the total cross sectional area of all the atoms in one clump while 
Q(27?uv)^ is the cross sectional area of one clump at r = 2/?uv- These can only be equal if clumps 
have the scale of individual atoms. Secondly, Eq. ( I30> can be expressed as 47r(27?uv)^ = Mcrs^ 
Here ^^ is the total area of all clump atoms per unit radial distance, so ^^^ is the total area of 
all clump atoms in a scale length 6. Thus, since 47r(2/?uv)^ is the spherical area at r = 2R\jv, the 
scale 6 defines the range of r around 2/?uv in a smooth wind over which the clump atoms just 
cover the sphere. 



4. Discussion and conclusions 

4.1. Conclusions regarding thick clumps 

Using a simple model we have shown that while clumping reduces the mass loss rate of WR stars 
required by radio emission measures, it also reduces the wind optical depth and hence multiple 
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scattering and momentum delivery. The nett result is that thick clumping worsens the momentum 
discrepancy rather than solving it. This is not the case for thin clumps, as we discuss below. 

4.2. Discussion of assumptions 

First we comment on various simplifications we have made which might modify our results 
somewhat. These include the approximation of constant u, relaxation of which does not seem 
very likely to change our results much since the radio emission measure is produced well out in 
the winds beyond where wind acceleration starts. Secondly, in common with many authors, we 
have so far assumed that the inner boundary /?,ad (and hence the volume) of the radio emission 
measure region does not change with clumping. In fact, one might expect, by analogy with the 
UV optical depth (Eq. ( I17» . that the radio optical depth might fall with increased clumping, 
reducing Ri-^d and increasing the radio source volume and emission measure. To check this, we 
first want to know the radio optical depth t' for a clumped wind. This is roughly given by 

r°° M 

T'(r)^ N(r)—o-'dr (31) 

J, m 

where cr' is the relevant cross section per proton. However, we have to note that the main radio 
absorption mechanism is free-free opacity which is density dependent (cr' cc n^ ~ Ijr^) and we 
have to write 

^' = cr;^^ (32) 



where n,. = rig, and cr' - cr'g are defined in any reference level r = r,,- Then Eq. OH becomes 



\2nv Ho \ml Q.Sr' ' 



.V)=^^(^)-^, (33) 

< n,, \m I 



so that now /?„(! is given by r'(/?rad) = 1, namely 



^r= 



rad 






Consequently, increasing M for given Q, S or decreasing Q, S for a given M (i.e., increasing 
dumpiness) actually makes /^rad bigger, not smaller, in a clumpy wind because of the density 
dependence of free-free absorption. 

The corresponding emission measure expression is now as before but based on the new 
clumping dependent value in Eq. ( I34t of 7?rad which leads to 

3'/^^ /^fV no \'^'/M\^/^ 



which does increase as we increase M/Q.6 (i.e., dumpiness) for a given M but now with e oc 
[MM/iQd)]-^^^, instead of [MM/(Q.6)] for the constant R^^d case. Thus although dumpiness still 
reduces M for a given e, it does so less than with constant R^^i and likewise, thick clumps are 
now even less helpful to the momentum problem. 
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4.3. Discussion of tliin clumps 

Throughout, we have considered here the case when each clump is an optically thick scatterer 
at line-driven wavelengths. The other case of relevance to reducing rj is where there is strong 
clumping (/ <k 1) but with very large numbers of very small scale clumps, each optically thin. 
We have considered here the case when each clump is an optically thick scatterer at line driving 
wavelengths but with a large overall wind optical depth. In this case our analysis does not really 
apply and the driving has to be described by multiple scatterings over the small Sobolev optical 
depths of many successive individual thin scatterers. The nett effect here is to yield the same 
momentum delivery rate as in a smooth wind but to reduce M for a given s by increasing the 
optically thin radio emission measure ei of each clump. In this thin clump limit, clumping does 
not reduce momentum dehvery by photon escape, as happens for thick clumps, but reduces M 
for a given e. As a solution to the momentum problem, this scenario puts the entire onus on 
the reduction of M/s and of 77 by a factor /'^^ and for the extreme case of 77 a; 100 requires 
/ X 10 '*. How such a huge compression of wind matter into clumps is achievable physically, 
and whether it can be done without making the clumps optically thick in the radio (bearing in 
mind that the radio opacity oc n^) are questions that must be addressed before clumping solutions 
of the momentum problem, are to be accepted. 

4.4. Otiier solutions to the momentum problem? 

Since clumping solutions, thick or thin, remain questionable, we briefly mention here some other 
possibly relevant factors. Some have advocated non-spherical models for the massive WR flows. 



Lamers & Pauldrach (I1991> developed a bi-stability model for early-type stars (e.g., B[e] stars), 
and Poe et al. ( I1989> proposed a two-component model with fast winds from polar regions and a 
slow flow from equatorial regions, later termed the Luminous Magnetic Rotator model (Cassinelli 



J1991» . In aspherical models such as these, the high radio flux arises from the denser equatorial 
region, whereas the high terminal speeds derive from a line-of-sight that lies perhaps in a broad 
polar region. These two-component structures have not been widely accepted because few WR 
stars show substantial (non-varying) intrinsic polarizations, as would be expected for stars that 
have a dense equatorial flow (Harries, Hillier, & Howarth ll998> . However, Taylor & Cassinelli 
(1992) studied the cancellation of polarization owing to a more tenuous polar wind, and it was 
surprisingly effective. Whether or not the solution to the momentum problem lies in the aspherical 
structure of such a rotational model, or yet some other aspherical picture, our paper has shown 
that the wide spread belief that the solution lies in the clumpy + multiple scattering picture is too 
simplistic. 
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Appendix A: Relation between Covering factor Y, number of scatterings A^„ 
and optical depth r 

Y is the fraction of the solid angle around a star that is covered by scatterers. Let A and to be 
the cross section and solid angle for one scatterer at r, so that w = A{r)/r^ and let Nsir) be the 
space density of scatterers, then the covering factor at r is the total solid angle of all the scatterers 
divided by An, namely 

F=— I N,(r)cj4-nr^dr^— I N Ar) -^ 4-n r^ dr ^ I N,(r)A(r)dr ^ T(r). (A.l) 

A photon travelling in a medium with typical size D, density Ns, and particle cross section A, 
on average undergoes A^, scatterings before escaping. Its mean free path / = IjiNsA), is related 
to Ns by I ^iNl = D, so 

A^.v = (j) ^{NADf^T^. (A.2) 
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